Universal three-body parameter in heteronuclear atomic systems 
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In Efimov physics, a three-body parameter (3BP), previously regarded as nonuniversal, uniquely 
defines bound and scattering properties of three particles. A universal 3BP, however, have been 
recently shown in experiments and theory in ultracold homonuclear gases. Our present study further 
predicts a universal 3BP for heteronuclear atomic systems near broad Feshbach resonances, and 
provides physical interpretations for its origin. We show that for a system composed of two light 
and one heavy atoms, the physical origin of the universal 3BP is similar to the homonuclear case 
while for systems composed by one light and two heavy atoms the universality of the 3BP is now 
mostly controlled by the heavy-heavy interatomic properties. This new universality is explained 
by the universal properties of the van der Waals interactions in a simple Born-Oppenheimer (BO) 
picture. Finally, we show the numerically determined 3BPs for some combinations of alkali atoms 
used in ultracold experiments. 

PACS numbers: 



The Efimov effect [1], predicted in the 70's by Vitaly 
Efimov, has attracted broad interest in studies in atomic 
and nuclear physics [2]. Motivated by the general in- 
terest to understand this esoteric quantum phenomenon 
and its consequences on ultracold quantum gases, signifi- 
cant progress has been made in studies of Efimov physics 
during the last decade. In particular, many theoretical 
predictions and new phenomena have been identified due 
to the extraordinary ability to control the interactions in 
ultracold quantum gases [3]. 

The Efimov effect is characterized by the formation 
of an infinity of trimer states (Efimov states) when the 
two-body s-wave scattering length a — > oo. The ener- 
gies of Efimov states, E n , follow a universal scaling law: 
E n = Eoe~ 2n7r ^ s ° , where sq is a universal constant that 
depends only on the mass ratios, the number of resonant 
pairs of interactions, and the identical particle symme- 
try. The ground Efimov state energy Eq, therefore, fully 
determines the Efimov spectrum and it is usually con- 
sidered as the three-body parameter (3BP). Other defi- 
nitions of 3BP exist, however, based on scattering prop- 
erties of the system. For instance, in ultracold quantum 
gases, the formation of an Efimov state can be mani- 
fested by the loss of atoms via three-body recombination 
via resonance (a < 0) or interference (a > 0) features in 
the loss rate K% as a function of a. The value of a where 
the first resonance or interference occurs, a*_ or ag, re- 
spectively, are also equally good 3BPs (relation between 
Eq, a*_ and can be found in Ref. [2]). Nevertheless, 
since the early years after Efimov's prediction, funda- 
mental assumptions led to the expectation that the 3BP 
should be nonuniversal [2]. 

Despite such expectations, ultracold experiments with 
alkali atoms [4-13] have observed a universal value 
for a*_. For homonuclear systems it was found that 
fll~- 9.1r v dw, where r v dw is the van der Waals radius 
[3]. The universal homonuclear 3BP was subsequently 
discussed by different theoretical studies [14-17]. In par- 



ticular, the adiabatic hyperspherical (hereafter "adia- 
batic") picture used by Wang, et al. [14] shows the exis- 
tence of an effective three-body repulsion when the hy- 
perradius R — the overall size of a three-body system 
— is about 2r v dw m homonuclear atomic systems, pre- 
venting three atoms from approaching shorter distances. 
The Efimov states are affected by this repulsion, yield- 
ing a value for the 3BP that is universally determined by 
r v aw [14]. Nevertheless, in a heteronuclear system the ex- 
tension of the universality in 3BP is not straightforward, 
particularly for the Efimov states with non-resonant in- 
teraction between one of the pairs. 

In this Letter, we predict universal 3BP in heteronu- 
clear atomic systems with two identical bosons (A) and 
one distinguishable atom (X). For a fixed mass ra- 
tio, the heteronuclear 3BP is found to depend only on 
r v dw between the pairs and the homonuclear scattering 
length. We, however, give different interpretations for 
the physical origin of the universality in the "Efimov- 
favored" systems (sq > 1) in the extreme of two heavy 
and one light atoms and the "Efimov-unfavored" sys- 
tems (so < 1) in the extreme of two light and one heavy 
atoms, respectively. The universality of the 3BP in the 
"Efimov-favored" systems is intuitively understood in the 
BO approximation, via the universal properties of the 
van der Waals interaction between the heavy particles 
at small distances. On the other hand, the universal- 
ity in the "Efimov-unfavored" systems will be shown to 
have similarities with the one found in homonuclear sys- 
tems [14]. Finally, the 3BPs for some experimentally 
available atomic species are listed in terms of a*_ and a^. 

The wavefunction for the relative motion of three 
atoms ^ is determined by the three-body Schrodingcr 
equation (in a.u.) 
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where r is the displacement vector between the two A 
atoms with mass M, and p is the vector from the center 
of mass of the A atoms to the X atom with mass m. 
Here we use the Lennard- Jones [18] potential that has a 
van der Waals tail to model the interactions between the 
atoms with distance taa (?ax)'- 
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The scattering length claa (a ax) is changed by tuning 
the short-range cut-off r C: AA (r c ,Ax)- We should note, 
however, that our single-channel treatment of the atomic 
interactions implies that the results presented here are 
valid for broad Feshbach resonances and less so for nar- 
row ones [19]. 

Equation (1) can be solved to a desired numerical accu- 
racy in the hypersphcrical representation [20, 21], where 
the hypcrradius R is defined by 



M 



2 Aim 



p 2 2M + m' 



(3) 



/i being the three-body reduced mass [22], and the set 
of hyperangles f2 defines the three-body geometry. Ex- 
panding VP in the complete, orthonormal adiabatic ba- 
sis [20, 21] $ v by 



* = J2F U!E (R)<Z> u (R;Q), 



(4) 



Eq. (1) reduces to a set of coupled hyperradial equa- 
tions: 
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with P vll i and Q vv > the nonadiabatic couplings. Typi- 
cally, for well-separated adiabatic potentials U V (R) the 
spectrum of the three-body system is determined by the 
effective potentials W V (R) = U v (R) — Q vl/ (R)/2^, and the 
corresponding channel wavefunctions F u .e{R)- 

For Efimov-favored cases (M/m ^> 1) the adiabatic po- 
tentials U y (R) are not isolated at small R, and the strong 
nonadiabatic couplings near the sharp avoided crossings 
make an analysis of 3BP based on W V {R) impractical. 
These crossings can be removed in the BO approxima- 
tion where r is taken as the adiabatic variable. 

Here we restrict our discussions only to the a g sym- 
metry [23] since the Efimov effect only involves s-wave 
interactions [1]. After a separable form of the wavefunc- 
tion is assumed, 



(0) 



the three-body problem reduces approximately to an ef- 
fective two-body problem, and the three-body spectrum 



can be determined by solving the BO channel wavefunc- 
tion (r) via 
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where U BO (r) is the BO potential, obtained by solving 
Eq. (1) for fixed values of r. It is well known that when 
|aytx| 3> f vc iw,AX, the potential relevant to the Efimov 
effect has the universal long-range behavior [2] : 
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with xo ~ 0.567143. In the BO limit, xo is connected to 
the Efimov scaling constant so by Sq rj XqM/2to — 1/4. 
The 3BP, however, is determined by the behavior of the 
potential U BO (r) near the short-range radius ro, where 
ro is the larger of r v aw,AX and r vt rw,AA- The numerical 
solutions of Eq. (1) shows that in this region U®°(r) ps 
Vaa(t) and since the 3BP is controlled by the position 
of the last node in (r) near r = ro [2] , the property 
of the potential Vaa(t) is the key to determine the 3BP. 
This is in strong contrast to the homonuclcar system [14] 
where the short-range details of the interactions are not 
referenced at all. 

From the universal van der Waals theory [24, 25] we 
know that the solution to the two-body Schrodingcr 
equation with the potential Vaa{i") is universally deter- 
mined by claa and r v dw,iA- The channel wavefunction 
F BO (r) for the Efimov states at short-range r < ro there- 
fore also shares these universal properties, as demon- 
strated in Fig. 1 for 174 Yb2 6 Li system with at unitarity 
(a ax = claa = oo). Here we assume a J n = + symme- 
try for the total orbital angular momentum J and total 
parity H Figure 1(a) shows that for fixed claa, decreas- 
ing the short-range cut-off r Cy AA (consequently creating 
more A2 bound states) only builds more oscillations in 
F®°(r) at small r without changing its long-range be- 
havior. In fact, the universality of the 3BP is supported 
by the position of the last node in F„j§ (r) (and therefore 
the corresponding energy), which changes only within a 
few percent when the depth of the potential VaaM is 
changed by more than 3 orders of magnitude for values 
°f r c,AA shown in Fig. 1(a). 

Although the BO approximation should yield the exact 
3BP in the limit M/m — > 00, it is important to know the 
significance of the non-BO correction in realistic atomic 
systems. The role of the A2 + X break-up channels in 
determining the 3BP is of particular interest due to the 
inability of the BO representation to describe such chan- 
nels. To this end, we solve Eq. (5) to obtain the exact 
Efimov spectrum. Figure 1(b) shows a comparison of 
the BO and the hyperspherical radial probability den- 
sities (|^°(r)| 2 and \F y (R)\ 2 ) for the first excited Efi- 
mov state at unitarity. The good agreement shown in 
Fig. 1(b) was observed numerically for mass ratio M/m 
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FIG. 1: (color online) The wavefunctions for the first ex- 
cited Efimov state at aAA~ciAx=oo. (a) The BO wave- 
function F?°(r) with increasing number of A2 bound states. 
The short-range cutoff is tuned at r c .AA/r v dw,AA ~0.92, 
0.57, 0.38, 0.31, 0.26 to give 1, 2, 4, 6, 8 A 2 s-wave bound 
states, respectively, (b) A comparison of the BO probabil- 
ity density |i ? I f°(r)| 2 and hyperspherical probability densities 
|i^(_R)| 2 for the first excited Efimov state at unitarity. Here 
we use A= 174 Yb and X= 6 Li, with r v dw,AA=78.7 a.u. and 
r v dw,Ax=38.1 a.u. 



ranging from 29 to 14, corresponding to Yb^Li and Rb^Li 
systems [26] . 

To study the non-BO correction to the Efimov energy 
spectrum more quantitatively, we first recall [24, 25] that 
the position of the last node in F®°(r) for a pure van 
der Waals interaction also depends on ciaa, suggesting 
a universal dependence of the 3BP on a a a- Figure 2 
shows such dependence for the low-lying Efimov states 
for the 174 Yb2 6 Li system. The spectrum is calculated 
by both the BO approximation and the hyperspherical 
representation. The energy of the ground Efimov state 
calculated by the BO approximation agrees almost per- 
fectly with the exact hyperspherical calculations. (The 
ground Efimov state is defined by the first Efimov state 
that appears at a ax = o*_.) The energies for the ex- 
cited states, however, start to deviate due to the finite 
mass ratio correction to the Efimov scaling factor sq in 
the BO calculations. Nevertheless, the overall agreement 
still provides further evidence for universality of the 3BP. 

The Efimov spectrum shown in Fig. 2 is periodic in 
1/a.AA, which implies that as the depth of the two-body 
potentials Vaa increases so that ciaa changes periodically 
from negative to positive, an Efimov state will follow the 
path in the spectrum and gets more deeply bound until 
it moves below the heteronuclcar van der Waals energy 
scale and becomes non-Efimov. 

The role of the A2 + X break-up channel is shown in 
Fig. 2 for the cases where ciaa ^S> Tq. Evidently, all the 
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FIG. 2: (color online) The Efimov energy spectrum E n (a,AA) 
for 174 Yb2 6 Li at the heteronuclear unitarity (clax =00). The 
"IBS" and "2BS" hyperspherical calculations have different 
r c ,AA and v C: ax to give 1 and 2 s-wave bound states for both 
A + A and A + X pairs, respectively. The error bar on the 
hyperspherical data points indicates the width of the corre- 
sponding Efimov state. The "Analytic" curve is calculated 
from Eq. (11) (see the discussion in the main text). 



excited Efimov states acquire large widths when the first 
excited state moves right above the A2 + X threshold. 
In fact, as will be discussed later, in general this thresh- 
old energy is a dividing point between two different re- 
gions of the Efimov energy spectrum: a lower Efimov 
spectrum controlled by larger Sq (corresponding to three 
resonant interactions) and a upper Efimov spectrum con- 
trolled by smaller so (corresponding to two resonant in- 
teractions) [27, 29]. But this division goes away in the 
BO limit because the values of sq and Sq becomes iden- 
tical [27]. 

The simple scaling behavior of the BO potential 
U®°(r) allows us to determine the 3BP for the Efimov- 
favored systems analytically. The key step to deter- 
mine the 3BP, in our case here Eq, is to find the po- 
sition of the first node in F®°(r) in the Efimov re- 
gion r > rrj. To this end, we divide the BO poten- 
tials in the different regions and use their approximate 
forms. The first region is the van der Waals region 
with r < r_ where Vaa(^) is dominant, and we let 
Uv°{r) = —Ce.AA/r 6 . The second region r_ < r < re- 
characterizes the deviation of the potential from both 
its short-range and asymptotic forms, and is approxi- 
mated as U®°(r) = —C e ,AA/r 6 — (cvxo) 2 /2mr 2 , where 
a 2 rj 2 from numerical observations. In the last region 
we let U®°(r) = — Xg/2mr 2 . The zero-energy solution in 
each region can be written down analytically in terms of 
Bessel's functions J v {x) and N v (x), and will be matched 
at the boundaries. For simplicity, we let r_ and ?*+ be 
nodal positions of these solutions, leading to the following 
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transcendental equations 
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As shown in Fig. 2, Eq. (11) gives a reasonable estimate 
of Eq when compared to the exact numerical results. 

For the Efimov-unfavored cases (M/m < 1) the BO 
picture becomes invalid. In contrast to the Efimov- 
favored cases, our hyperspherical calculations now show 
strong nonadiabatic corrections, suggesting a different 
physical regime for the universality in the 3BP. The en- 
ergy spectrum is now well separated into two parts by 
the A 2 + X threshold energy when \clax | ^> a-AA ro- 
The upper and lower parts of the spectrum follow the 
Efimov geometric scaling for two- (sq) or three- resonant- 
interaction (sq > 50)7 respectively [27, 29]. Therefore in 
contrast to the Efimov-favored cases there are two 3BP to 
be considered. The two distinct 3BP are manifested by 
the two well-separated Efimov potentials W U (R) which 
become more weakly coupled as M/m decreases. Specif- 
ically, the upper Efimov spectrum comes from the three- 
body potential associated to the AX+A channel, whereas 
the lower spectrum originates from the potential associ- 
ated to the A2 + X channel. 

To simplify our discussion, the analysis below is made 
only for a ax = 00, however, same conclusions hold for 
finite values of a ax- Figure 3(a) shows the universal 
behavior of the AX + A Efimov potential W V (R) for 
133 Cs2 87 Rb. The universal form of the potentials after 
a scaling by claa strongly suggests a universal 3BP that 
depends only on a a a- Figure 3(b) shows the universal 
properties of the A2+X Efimov potential, which resemble 
those in the homonuclear case [14]. In this case, W V (R) 
also displays a repulsive barrier near R > ro independent 
of the number of two-body bound states . It is also impor- 
tant to note that when ciaa <C ro there is only one contin- 
uous Efimov spectrum with scaling constant Sq and the 
Efimov potential W U (R) shares the same universal fea- 
tures as the potentials shown in Fig. 3(b). Although the 
Efimov potentials are universal for the Efimov-unfavored 
cases, we should point out that the exact shape of the po- 
tentials near their minimum depends on the mass ratio 
M/m. 

Table I summarizes our results obtained for the 3BPs 
relevant to the three-body recombination process, A + 
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FIG. 3: (color online) The adiabatic hyperspherical poten- 
tials W v (R) for the AX + A channel (a) and the A 2 + X 
channel (b) for an Efimov-unfavored system with A= 133 Cs, 
A'= 87 Rb. Here r v dw,AA = 101 a.u. and r v dw,Ax=90.9 a.u.. 
The aAA labels in the plots are also in a.u.. 



so 



ciAA.bg (a.u.) ag (a.u.) a*L (a.u.) 



174 Yb 2 6 Li 
133 Cs 2 6 Li 
87 Rb 2 6 Li 
41 K 2 6 Li 
23 Na 2 6 Li 
87 Rb 2 40 K 

133 Cs 2 87 Rb 0.535 1.060 
41 K 2 87 Rb 0.246 0.961 



2.246 2.382 104 [32, 33] 
2000 [34] 
100 [35] 
62 [36] 
100 [37] 
100 
2000 
62 



1.983 2.155 
1.633 1.860 
1.154 1.477 
0.875 1.269 
0.653 1.125 



1.3 x 10 d 

9.6 x 10 2 
3.8 x 10 2 

3.7 x 10 2 
1.5 x 10 3 



-8.4 x 10 2 
-1.4 x 10 3 
-1.6 x 10 3 
-2.4 x 10 3 
-1.3 x 10 4 



2.8 x 10 3 < -3 x 10 4 
2.3 x 10 3 < -4 x 10 4 
> 7 x 10 3 < -1 x 10 6 



TABLE I: The universal Efimov scaling constants so, Sq and 
the 3BPs clax = ag and ciax = aL obtained by keeping aAA 
fixed at its background value (a,AA,b g )- 



A + X -> AX + A [27, 28], for some combinations of 
alkali atoms. The 3BPs are given in the form of the ex- 
perimental observables al and Oq in a ax near a A + X 
Feshbach resonance. It is important to note that the 
knowledge of these universal 3BPs permits more precise 
calibration of the position of a Feshbach resonance in ul- 
tracold quantum gases, providing great benefits to more 
accurate measurements on strongly interacting quantum 
gases. Finally, we have also tested the universality of 
the 3BPs against three-body forces [30, 31] and have 
found a negligible effect on the 3BPs. It has been pre- 
dicted [28] that al/aQ = e 7T ^ 2s ° in the zero-range limit 
with aAA = 0. However, due to the complicated inter- 
play of those finite length scales in realistic systems [26] , 
significant deviations of this relation are observed in the 
numerical results in Table I. 

It is also worth mentioning that our results for Rb2K 
and estimations on 3BP for K^Rb do not agree with 
the experimental results by G. Barontini et al. [38]. We 
believe such disagreement might be related to effective- 
range corrections and requires further investigations. 
Our Rb2K results, however, are more consistent with 
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JILA experiment [39], where no Efimov resonances were 
observed for a ax > — 10 4 , although a final conclusion 
cannot be made due to the thermal saturation in that 
experiment. 

To summarize, we have predicted and analyzed univer- 
sality in the 3BP in heteronuclear atomic systems. The 
physical origin of this universality, however, is different 
depending on the mass ratio between the atoms. For the 
Efimov-favored cases, the universality is explained by the 
universal property of the atomic van der Waals interac- 
tion between the heavy atoms, which yields an analyti- 
cal determination of the 3BP in the BO limit. For the 
Efimov-unfavored cases, the universality is manifested by 
a universal effective three-body repulsion in the range of 
the van der Waals radius shielding three atoms from the 
complicated interactions at shorter range. Finally, the 
3BPs we have calculated can be used in ultracold exper- 
iments for calibration of the positions of the Feshbach 
resonances. 
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